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Abstract. Simulated color-magnitude diagrams are used to investigate the recent star formation history in 
NGC 5128. The comparison of the simulations with the observed UV color-magnitude diagram for a field in 
the north-eastern shell, where recent star formation is present, constrains the initial mass function (IMF) and the 
duration of the star formation episode. The star formation burst is still on-going or at most has stopped some 2 
Myr ago. The look-back time on the main sequence is set by the incompleteness of the f/-band observations at 
V ~ 26, and is of the order of 50 Myr. The post main sequence phases have a look-back time up to 100 Myr, 
setting the maximum observable time for the star formation in this field. The comparison of the observed and 
simulated luminosity functions and the number of post main sequence vs. total number of stars favours models 
with active star formation in the last 100 Myr. The data are best fitted by an IMF with Salpeter slope (a — 2.35), 
and are also consistent with slopes in the range of 2 ^ a ^ 2.6. They exclude steeper IMF slopes. The mean 
star formation rate for a Salpeter IMF in the range of masses 0.6 < M < 100 M© within the last 100 Myr is 
1.6 X lO"''M0yr-^kpc"^ 

Key words. Galaxies: elliptical and lenticular, cD - Galaxies: stellar content - Stars: fundamental parameters - 
Galaxies: individual: NGC 5128 

1. Introduction served in these galaxies is an important complementary 

, ; approach to the studies of the high redshift Universe, aim- 

. The two main traditional models of elliptical galaxy for- j^g to address the questions of elliptical galaxy formation 
■ ■ ' mation, the "monolithic collapse" (Partridge & Peebles g^^j evolution 
I1967|l and the "hierarchical assembly" model (Toomre 
I1977|l assume very different mechanisms and epochs of 
spheroid assembly. Correspondingly, they also predict dif- 
ferent mean ages and metallicities for the stars that build 
these galaxies. Until the advent of the HST and 10m class 

ground-based telescopes the studies of the stellar popula- At the distance of 3.8 Mpc (Soria et al. 119961 Rejkuba 
tions in the elliptical galaxies were restricted to integrated |2004f) NGC 5128 (=Centaurus A) is the dominant galaxy 
light photometry and spectroscopy. In the last few years, of the nearby Centaurus group and one of the closest giant 
however, the nearest E/ SO gal axies have been resolved into elliptical galaxies to us. It is also one of the nearest AGNs 
stars (e.g. Soria et al. |1996' Davidge & van den Bergh and radio galaxies. It presents a number of pecuharities: 
120011 Schulte-Ladbeck et al. 2003 ). The study of the star the central dust lane with HII regions and OB associations 
formation histories of the resolved stellar populations ob- at its edge, and the stellar shells which have presumably 

been formed in a recent accretion of a smaller companion 



o 
o 

(N 

in 



> 

o 

m 

o 
o 

O 



1.1. Recent Star Formation in tine halo of NGC 5128 



Send offprint requests to: M. Rejkuba galaxy (see review by Israel Moreover, young blue 

* Based on observations collected at the European Southern supergiants have been identified associated with the opti- 

Observatory, Paranal, Chile, within the Observing Programme cal fi laments i n the north-eastern halo (Graham & Price 

63.N-0229 [TOHTl Graham [TUMji . 
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The young stars and OB associations in the halo of 
NGC 5128 are roughly aligned with the ionized gas and 
the radio jet (Morganti et al. [TMl Mould et al. 1^00171 
Fassett & Graham iZPOtil Rejkuba et al. 120011 12001 over 
several kiloparsecs. The preferred mechanism for this re- 
cent star formation in the halo at > 14 kpc away from 
the nucleus is through the direct interaction of the radio 
jet with the gas. The necessary amounts of neutral gas 
for this star formation have been detected in nearby fields 
through HI (Schiminovich et al. 119941 and CO molecular 
transitions (Charmandaris et al. pOOO). Charmandaris et 
al. suggested that the molecular clouds in the distant halo 
are associated with stellar shells and are thus remains of 
the last accretion episode in which NGC 5128 has merged 
with a gas rich companion galaxy. While the gas could 
have just "hung around" in the galaxy halo, the presence 
of the radio jet triggered the star formation in highly col- 
limated areas. The investigation of the star formation his- 
tory of the young stars and the determination of the age 
spread in the young stellar population can constrain the 
formation of the radio jet or the time of the last accretion 
event. The metallicity of the gas of the accreted galaxy 
should be reflected in the newly formed stars. 

The comparison of the optical color-magnitude dia- 
grams (CMDs) with stellar evolutionary isochrones and 
young LMC clusters has indicated that these young stars 
can be as young as 10-15 Myr old (Mould et al. 1^0001 
Fassett & Graham 151100!, Rejkuba et al. I2001|l . However, 
except for a comparison with one synthetic CMD (Mould 
et al. I2()()()|l , no detailed studies of the star formation his- 
tory in the halo of NGC 5128 were made. 

1.2. The Method 

Containing stars that were born during the whole life- 
time of the galaxy, CMDs represent snapshots of their 
star formation histories. Even when limited to the bright 
stars, the CMDs contain useful age and metallicity indi- 
cators: the colors of the red giant branch (RGB) stars are 
very sensitive to metallicity; the presence of bright asymp- 
totic giant branch (AGB) stars signals an intermediate-age 
population; the colors and magnitudes of the bright main 
sequence (MS) stars and blue supergiants depend on the 
most recent star formation history, as well as on the metal- 
licity. 

The most effective way to analyse the distribution of 
stars in the CMD and to quantify the star formation his- 
tory of a galaxy is through the comparison of the observed 
CMD with synthetic CMDs computed using stellar evolu- 
tion models. This technique is based on the simulations 
of a large number of CMDs constructed assuming differ- 
ent star formation rates (SFR(i)), chemical enrichment 
laws {Z{t)) and initial mass functions (IMF). The photo- 
metric errors and the incompleteness of the data are also 
simulated. Each simulated CMD then depends on all the 
main parameters which determine the distribution of stars 
in a CMD. In this way it is best accounted for the evo- 



lutionary behaviour of stars, small number of stars, and 
spread due to observational effects. The synthetic CMDs 
are then compared to the observed ones in order to de- 
termine the acceptable range of parameters and the best 
fitted models. This technique is currently used by several 
groups to investigate the star formation history (SFH) in 
dwarf galaxies (Ferraro et al. ll989l Tosi et al. 119911 Bertelli 
et al. T olstoy IT!?^ Galla rt et al. Il99f)al I1996bl 

DolphinTW, Greggio et aLTOM, Hurley-Keller, Mateo & 
Nemec il998i Hernandez, Valls-Gabaud & Gilmore ITOMl 
Cole et alHOOOl Olsen lT^ Schulte-Ladbeck et al. 120001 
Harris & ZaritskyEOOB- 

Here we apply this technique to investigate the recent 
star formation history of the halo of NGC 5128. The ob- 
servations were taken with VLT using FORSl and ISAAC 
in the U, V and Ks bands. Their full description as well 
as data reduction and the resulting CMDs were presented 
in Rejkuba et al. ip]UT|l . 

This paper starts with the description of the synthetic 
CMD simulator, including a detailed description of the al- 
gorithms used. In the next section the young stellar popu- 
lation in the north-eastern halo field is simulated. We test 
two different metallicities and a range of star formation 
histories. The initial mass function is also investigated. 



2. Description of the synthetic CIVID simulator 

Our synthetic CMD simulator is developed based on the 
code used by Greggio et al. I|1998|l . The construction of a 
synthetic CMD goes through the following steps: 

1. Random extraction of mass according to the initial 
mass function (IMF), age according to a chosen star 
formation rate (SFR) law, and metallicity according to 
the chemical enrichment law [Z{t)). 

2. Positioning of the extracted star as a function of mass, 
age, and metallicity on the theoretical HR diagram 
(log(L/L0), log(Tcff)) interpolating within stellar evo- 
lutionary tracks. 

3. Transformation of the theoretical position of the star 
in the HR diagram (log(L/L0), log(Teff)) into the ob- 
servational CMD assigning the bolometric corrections 
through interpolation within conversion tables. 

4. Application of observational uncertainties, photomet- 
ric errors and incompleteness, according to the results 
of artificial star experiments on the real data. 

5. Iteration until a certain criterion (e.g. the total number 
of stars generated, the total mass extracted, or the to- 
tal number of stars in a certain region of the simulated 
CMD) is satisfied. 

The algorithms used for steps 1 and 2 are relatively 
straightforward, and used in a similar way by most groups. 
Hence we describe them in the Appendix. In the following 
sections we address steps 3 and 4 in more detail as the 
final results depend critically on them. 
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Fig. 1. Left: Fagotto et al. 1] 1994(1 tracks for Z = 0.004, transformed into the V vs. U — V plane by applying bolometric 
corrections obtained by convolving spectral energy distributions of Kurucz stellar atmospheres (see text) with the two 
sets of U BV filters: Landolt (1 19921 red lines with crosses) and Buser & Kurucz 1)19781 blue lines with triangles) The 
masses of the tracks are indicated. Right: The same as on the left, but for / vs. V — I. The V filter transmission curves 
are from Landolt H1992I red lines with crosses) and Buser & Kurucz p978l blue lines with triangles). The / passband 
is the same in both cases. 



2.1. Bolometric corrections 

Bolometric corrections (BCs) are needed to convert the 
theoretical log(L) — log(Tcfr) point into magnitude and 
color, to be compared with the data. An accurate knowl- 
edge of the photometric passbands is as important to syn- 
thetic photometry as the knowledge of the standard stars 
and the proper calibration procedure is to observational 
photometry. The calibration of our data was described in 

mm 



detail in Rejkuba et al. 1)2001(1 . We have derived a full set 
of calibration equations including the color term for the 
optical photometry, while only the zero point was derived 
for the near-IR photometry. Thus our photometric system 
consists of Landolt ijTUU^ U and V and ISAAC Ks (Cuby 
et al. 12003). 

The Kurucz ijTUMjl database of ATLAS9 synthetic 
spectra covers a large range in temperatures (3500 < 
Teff < 50,000 K), and gravities (0 < \ogg < 5). Recently 
these models have been re-calculated by Castelli et al. 
1)1997)1 and carefully compared to stellar data. Of the two 
sets in Castelli's grid, the one without overshooting is 
found to yield the best match to the observations (Bessell 
et al. 1998). Therefore we adopt the "NOVER" files 
at |http:/ /kurucz. harvard.edu/ for metallicities [M/H] = 
—2.0, —1.5, —1.0, —0.5, 0.0 and-|-0.5, and microturbulent 
velocity ^ = 2 km s~^ to derive BCs as function of gravity, 
effective temperature, and metallicity. This is done by con- 
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Fig. 2. Comparison of the Landolt 1)19921 full lines) and 
Buser & Kurucz 1)19781 dotted lines) UBV passbands. 



volving the spectral energy distributions from the Kurucz 
stellar library with filter transmission curves S\: 







dX 



^ [http:/ /www.eso.org/instruments/isaac/index.html| 
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etry in the V and i^j-bands. In the UV only bright young 
stars are observed, the giants are too faint in the J7-band. 

We calculate the log 5 of each simulated star with 
(log M/M0, log L, Toff) using the following two equations: 



4 000 5000 
X 



Fig. 3. The spectral energy distributions of the Sun and 
Vega are compared with the Landolt ( 199^) U BV pass- 
bands (full lines) and Buser & Kurucz (,1978|il U filter (dot- 
ted line). 



The magnitude and the bolometric correction are then by 
definition: 

M, = -2.51og(F,) + ZP(2) 
BQ = Mboi - M, = -2.5\og{aT^ff) + 2.51og(F,) - ZP{3) 



where aT^^j ^ j 



The calibration of the zero points in each band follows 
from the definition of the 1^-band bolometric correction 

—0.07, and the colors of Vega 



of the Sun to be BC^ 



M 



V ega 



M 



V ega 



0. The spectral energy distributions of 



the Sun and Vega were taken from the Kurucz library as 
well. 

It is known that the Kurucz H1993II spectral energy 
distributions are not very accurate for cool stars due to 
inadequate molecular opacities, non-LTE, and spheric- 
ity related effects (e.g. Morossi et al. I1998|l . To extend 
our bolometric corrections to low temperatures (Toff ^ 
4200K) , we used the empirically determined temperatures 
and bolometric corrections from Montegriffo et al. H1998(l . 
Montegriffo et al. do not have fully developed metallicity 
classes, but rather define metal-poor ([M/H] < —1.0 dex) 
and metal- rich ([M/H] > —1.0 dex) samples. The differ- 
ences in bolometric corrections between the two groups 
span a wide range 0.05-0.4 mag/dex. We also notice that 
their bolometric corrections are determined only for giant 
stars. However, the last issue is not a problem, since the 
low temperature stars in our case are giants with photom- 



L 



log(5/<?0) = log(M/M0) - 21og(i?/i?G 



(5) 
(6) 



where logg© = 4.437, ct = 5.6710-5erg/cmVK4, Lq = 
3.8261O33erg/sandR0 = G.gSggiO^^cm. The linearly in- 
terpolated bolometric correction, as a function of logTcft 
and log (7, is then applied to obtain magnitudes for the 
simulated stars. 

In Fig.^we compare the evolutionary tracks of Fagotto 
et al. l(TnM|l for Z = 0.004 transformed using two different 
sets of UBV filters. In the right panel the comparison of 
I vs. V — I shows an excellent agreement of the colors of 
tracks for both sets of filters, but in the left panel it is 
obvious that our bolometric corrections calculated using 
the Landolt filters produce wider range of colors for the 
tracks in U — V. We investigated in detail the origin of 
this difference, finding that it is entirely due to different 
U filter used in calculating the bolometric corrections. 

The Padova group uses Buser & Kurucz (|1978|l U BV 
transmission curves, while our photometry has been cali- 
brated to the Landolt (,1992;) system. The comparison of 
the passbands of the two photometric systems is given in 
Fig. El Note the large difference in the L/-passband wave- 
length coverage. The fj-band is known to be the most dif- 
ficult to reproduce due to the presence of the atmospheric 
cut-off, the presence of the Balmer break within that band, 
and the strong blanketing by metals (see Fig. The 
differences between the B and y-band filters are much 
smaller. Buser & Kurucz 1)1978(1 UBV system has two B 
filters, B2 and B3. Both of these systems claim to be based 
on the Johnson & Morgan H1953|l standard system with 
Buser & Kurucz U — B and B — V colors defined as follows: 



{U - B),yn = {U- B2) - 1.093 
[B ~ V),y^ = (B3 - y) + 0.710 



(7) 
(8) 



We have calculated the bolometric corrections by con- 
volving the Kurucz stellar library with the UB2B3V filters 
of Buser & Kurucz H1978|l . There is an excellent agreement 
between the transformed tracks and Padova isochrones 
(Bertelli et al. 1994J if the following definition of the U — V 
color is adopted: 



{U - V)syn = {U- B2) + {Bs - V) 



(9) 



The difference in the zero point is due to different cali- 
bration of the zero point bolometric corrections adopted 
by Buser & Kurucz with respect to the Padova isochrones 
(Girardi et al. 12002)1 . The Padova group adopted the stan- 
dard calibration of the synthetic photometry which is 
based on the colors and magnitude of Vega, similar to 
our procedure. 

This example shows how important it is to know pre- 
cisely the transmission curves and the characteristics of 
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the photometric system used. The cahbration of the syn- 
thetic magnitude system has to be done in the same way 
as the data, by convolving the spectral energy distribu- 
tions from the stellar library with the same filters used for 
the real observations. 

Our data have been calibrated on Landolt's stellar cat- 
alogue (Landolt n'992(l . but observed with the FORSl tele- 
scope+CCD-l-filter system, so that our system's through- 
put in the U band is indeed much closer to the Bessel 
U-filter (with a passband similar to the Buser & Kurucz 
I1978|l . mainly due to the CCD response curve. The best 
approach to calibrate our photometry would have been 
to use a standard system which matches the one used in 
the observations fBessell fl995|l . Unfortunately, our list of 
observed standards has no star in common with Bessell's 
l(T^ hst. 

The JJ-magnitude difference between the two systems 
depends on the star's spectral energy distribution, and 
therefore on temperature, gravity and metallicity. To a 
first approximation, the calibration equation corrects for 
this systematics as a function of the {U — V) color. 
However, a residual systematic error exists, which stems 
from applying to all of the observed stars the calibration 
regression valid for the standard stars. In order to estimate 
the residual error we have calculated the f7-magnitude dif- 
ference in the two U passbands of the complete grid of 
Kurucz models. The (FORSl - Landolt) difference of the 
U band Bolometric Corrections increases from ~ —0.2 to 
~ -1-0.4, going from the hottest to the coolest models in 
the grid. Once the color equation applied to the data is 
taken into account, the residual difference ranges from ^ 
-0.06 mag for the blue stars {U — V < 0), to ^ 0.25 mag at 
{U — V) > 2. We notice that this is the maximum residual 
error that we expect in our photometry, as it includes a 
very wide range of stellar parameters. Nevertheless, our 
photometry will tend to be too faint for the blue stars 
and too bright for the red stars. We will discuss in the 
last section the impact of this unceratinty on our results. 

2.2. Observational biases 

The synthetic CMD needs to be "corrected" for observa- 
tional biases before comparing it with the observed CMD. 
There are two main types of observational biases: 

— photometric errors: the magnitudes of stars are mea- 
sured with some uncertainty 

— incompleteness: some stars are lost due to crowding or 
high background 

Blending is an additional bias that is important in crowded 
fields, and that combines properties of photometric errors 
and incompleteness. 

We used the extensive artificial stars experiments (the 
full description can be found in Rejkuba et al. I2001|l to 
assign the observational biases. In the artificial star ex- 
periments the input magnitude is given, and the output 
magnitude measured at each position of the input mag- 
nitude. If the star was lost in the background noise, or 



Table 1. Fraction of stars in different parts of observed 
CMD - all stars, blue {U — V < 0), intermediate color 
{0 < U -V < 1.5), and red stars {U - V > 1.5) - that 
survive the error selection at maximum error of 0.1, 0.2, 
and 0.3, with respect to selection with a < 0.5. Only the 
stars with C/-band completeness limits larger than 50% 
are selected, corresponding to y-band magnitudes smaller 
than 25.5 in the bluest part of the CMD 



a cut 


all stars 


blue 


intermediate 


red 


a < 0.3 


81 % 


97% 


70% 


75% 


a < 0.2 


66 % 


87% 


53% 


55% 


a < 0.1 


37 % 


52% 


29% 


26% 





Fig. 4. Color-magnitude diagrams of the young stellar 
populations in the shell field of NGC 5128 with differ- 
ent cuts applied on the size of the photometric errors in 
both filters U and V. 



if it was completely blended with a much brighter object 
producing a magnitude difference between the input and 
the measured magnitude of more than 0.75, we considered 
it to be non-detected. Observational biases are applied to 
magnitudes to which the appropriate distance modulus 
and extinction were added. Thus, when applying the ob- 
servational biases we pick a random object from the input 
table of the artificial star experiments that has magnitude 
and color similar to the simulated star. The simulated star 
then receives a realistic photometric error (magnitude dif- 
ference between the recovered and input magnitude in the 
artificial star experiment), and incompleteness (the star 
in the artificial star experiment was detected or not). 



3. Young stellar population observed in the UV 
CMD 

The observations of the young stars in the north-eastern 
shell field of NGC 5128 consist of pairs of U and V-band 
images taken at ESO Paranal observatory UTl Very Large 
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Table 2. Main sequence and total life-times of the stars 
with masses larger than 4 Mq and metallicities Z — 0.004 
and Z = 0.008 (Fagotto et al. iraMji . 



M(Mg) TMS (Myr) rt,i (Myr) tms (Myr) rt<,i (Myr) 
Z=0.004 Z=0.008 



5.0 


101.89 


112.09 


104.91 


116.05 


6.0 


70.98 


76.74 


71.63 


77.71 


7.0 


52.58 


56.39 


52.40 


56.26 


9.0 


33.07 


35.25 


32.42 


34.52 


12.0 


20.20 


21.48 


19.85 


21.08 


15.0 


14.35 


15.27 


14.03 


14.93 


20.0 


9.805 


10.48 


9.44 


10.10 


30.0 


6.37 


6.86 


6.18 


6.66 


40.0 


5.01 


5.44 


4.89 


5.31 


60.0 


3.85 


4.20 


3.84 


4.23 


100.0 


3.11 


3.38 


3.12 


3.44 


120.0 


2.98 


3.24 


3.01 


3.33 



Telescope with FORSl instrument (Szeifert ElEF) The 
photometry is 50% complete for magnitudes brighter than 
V = 27.5 and U = 25 (Rejkuba et aL [M7T|l . The shal- 
lower [/-band photometry drives the incompleteness and 
the photometric errors in the UV CMD (Fig. ^ and thus 
limits our analysis to magnitudes brighter than V ~ 25.5. 

In Fig. m we display the V vs. U-V CMDs of this field 
with different values of DAOPHOT errors in both U and V 
applied as a selection criterion. The vertical blue sequence 
at [/ — ~ — 1 is populated with massive MS stars in 
NGC 5128, while most of the red stars belong to the fore- 
ground Galactic population. The diagonal cut-off at faint 
magnitudes is a consequence of the J7-band incomplete- 
ness and the sharp bright magnitude cut-off at F = 20 
is due to the saturation limit of the F-band images. The 
CMD with all stars that had a < 0.5 is presented in the 
upper left panel. Comparing it with the other panels it is 
evident that many of the red stars have large photometric 
errors. This is clearer in Tab. ^ where we list fraction of 
stars that remain in the CMD when a cut is applied at a of 
0.3, 0.2 and 0.1. Even when only the stars with the small- 
est photometric errors are selected, more than 50% of the 
blue stars survive and the young blue sequence remains 
well populated. We shall use this feature to investigate 
the recent star formation history in this field. 

Fig. O shows the observed CMD with the transformed 
Padova tracks superimposed, for Z = 0.004 (left) and 
Z — 0.008 (right). No spectroscopic estimate of the metal- 
licity of the young stars in Cen A is available in the lit- 
erature. The two considered Z values are motivated by 
the Rejkuba et al. H2001fl estimate, which is based on the 
best fitting isochrones. Notice, however, that due to the 
discussed difference between the Landolt and the Buser 
U filter, the Rejkuba et al. estimate needs to be recon- 
sidered with the appropriately transformed tracks. A dis- 
tance modulus of 27.8 mag (Soria et al. ll996|l in combina- 
tion with a color excess of E(J7— F) = 0.62 yields a good fit 

|http:/ /www. eso.org/instruments/fors/userman/, 



to the position of the MS stars in Fig.jS] According to the 
Cardelh et al. ifTM^ extinction curve, E{U - V) ^ 0.62 
implies E(B — V) — 0.35, which is much larger than the 
foreground extinction ^{B — V) = 0.11 from Schlegel et al. 
(1998) map, and implies a substantial amount of internal 
(perhaps differential) reddening. Rejkuba et al. estimated 
E{B -V) = 0.15, due to the different U filter. Choosing a 
somewhat larger distance modulus of 27.92 mag (Harris et 
alUnnni Rejkuba IMHIi does not change our conclusions 
significantly. 

The magnitudes and colors of the stars belonging to 
the Milky Way, that are expected to be found in a 6.'8 x 6'.8 
field of view of FORSl, at I = 309.5°, b = 19.5°, have been 
simulated using Besangon galactic model of stellar popu- 
lations (Robin & CrezeT^HH", Robin et al.]TM^- Realistic 
photometric errors and incompleteness has been applied 
through crowding simulations. In total, at least 340 stars 
with colors U — V > are expected to belong to the fore- 
ground Milky Way population (see Rejkuba et al. I20()l|l . 
They are overplotted over the observed CMD in Fig. [3 
as (green) filled triangles. Only a small number of stars 
with {U — V) > are expected to be blue and red su- 
pergiants evolving along the blue loops in the NGC 5128 
halo. Unfortunately, the width of the blue loops cannot be 
used to constrain the metallicity both because of the huge 
amount of the foreground contamination and, because of 
the uncertain extension of the blue loops fRitossa ll996|l . 

The location of the post main sequence (PMS) stars, 
blue and red supergiants, and the extension of the loops 
depends critically on the convection criteria or the mix- 
ing efficiency, the metallicity, the opacities at intermediate 
temperatures, the mass loss and the overshooting (Renzini 
et al. 119921 , Ritossa ri996|l . At larger metallicities the ther- 
mal conductivity of the envelope is lower leading to a 
larger thermal imbalance of the stellar envelope and to 
a faster expansion to the red. The more metal-rich models 
have thus redder red giant minima. The models that adopt 
efficient mixing in the semiconvective region during the 
hydrogen shell burning spend most of the PMS life-times 
as blue supergiants, while the less efficient mixing drives 
the runaway expansion toward the red fRitossa [TUM|l . The 
convective core overshoot during the MS in combination 
with efficient mixing can, however, produce models with 
wider loops (e.g. Fagotto et al. 11994(1 . 

The overshooting and opacity, constraining Tog, drive 
the amplitude of blue loops and, through the depen- 
dence of bolomctric corrections on the effective temper- 
ature (BC(Tcff)), they produce the sharp edge in the ob- 
served CMD. In particular the 9 Mq model for Z = 0.004 
has a very short loop (see Fig. |SJ . The blue excursion of 
this model is terminated at low temperatures, producing 
a sharp "red hook" in the simulated CMD (see below). 

The extremely short evolutionary life-time for stars be- 
tween the MS turn-off and blue edge leaves a "hole" in the 
CMD. On the other hand this area (-0.5 <, U-V 
and V ^ 23) of the observed CMD is populated. The pos- 
sible reasons include different metallicity, inadequate in- 
put physics (combination of overshooting and opacity pa- 
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Fig. 5. Comparison of the the UV CMD with stellar evolutionary tracks for Z = 0.004 (left) and Z = 0.008 (right 
panel) and masses of 6, 9, 12, 20 and 40 Mq. Overplotted are also the expected foreground stars (green filled triangles) 
from the Besangon simulations (see Rejkuba et al. l2U01|l . Along the line connecting MS turn-offs (large cyan dots) and 
along the line connecting blue edges (large black squares) most of the stars are expected according to the models. A 
box on the main sequence between magnitudes 23 and 24 indicates the area used to constrain the simulations. 



rameters), differential reddening, bolometric corrections, 
photometric errors and blends of several lower mass stars. 
Thus we expect not to be able to reproduce well this zone. 

The evolutionary lifetimes of the plotted tracks are 
listed in Tab. 12 Comparing the tracks to the data we 
see that the oldest MS stars in our observed CMD are 
~ 40 - 50 Myr old for the Z = 0.004 and - 50 - 65 Myr 
for the Z — 0.008 models. The faintest evolved stars are 
somewhat older (^ 80 Myr), but buried in the foreground 
contamination. We can thus derive the most recent star 
formation, and check the IMF for stellar masses larger 
than 7 M0 from the bluest stars in the CMD. 

Our simulations are constructed under the requirement 
of reproducing the stellar counts in a box on the MS, which 
contains 57 stars with 23 < V < 24 and U ~V < -0.7 
(Fig. EJ. The faintest stars in this box are 15 Myr old. 
Therefore, any SFH with the last star formation episode 



ceasing earlier than 15 Myr ago would not populate this 
box, and we do not consider these models. Given the short 
age range sampled by the blue stars, and given their rela- 
tively smooth distribution, we consider only episodes with 
constant SFR, and vary the starting {Tmax) and ending 
{Tmin) epochs of activc SF. We compare the synthetic and 
observed CMDs directly as well as the luminosity function 
along the MS for the stars brighter than V = 26, which 
are not affected critically by photometric errors and in- 
completeness. In the next sections we address the follow- 
ing questions: "How old are the oldest and the youngest 
stars in the observed UV CMD?" and "What is the slope 
of the initial mass function?" 
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Fig. 6. Synthetic CMDs for Z = 0.004 using the Salpeter IMF (a = 2.35), the appropriate distance, and a constant 
SFR in the interval that is indicated in each panel (beginning-end of the burst). For comparison the observed CMD is 
also displayed. In the lower part simulated luminosity functions (solid lines) are compared with the observed luminosity 
functions (dotted lines). In the observed and simulated luminosity functions the color and magnitude cuts are applied, 
and only stars bluer than (IJ — V) — —0.7 and brighter than V — 26 are compared. 



3.1. Age of the young stellar population in the 
NGC 5128 halo field 



The range of models probed is listed in Tab. El In the first 
column there is the model name, and in columns 2, 3 and 4 
the input parameters regarding the IMF and the duration 
of the star formation are given. The output total extracted 
mass in the range of 0.6 to 100 Mq, the mean SFR for 
the same mass range, and the explored interval of the star 
formation (At), as well as the total number of stars above 
the 80% incompleteness limit in F-band {V < 25.5), and 
the number of simulated blue plume stars {U — V < —0.7 
for V < 25.5) are listed in columns 6, 7, 8 and 9 for Z = 
0.004, and 10, 11, 12 and 13 for models using evolutionary 
tracks for metallicity Z ~ 0.008, respectively. Due to the 
small range of ages and to the smooth appearance of the 



MS, a constant star formation was considered. Therefore 
the columns 6 and 10 list: 



(SFR) = 



M, 



A7 



(10) 



Four different IMF slopes have been used in the simula- 
tions. The dependence of the simulations on the IMF slope 
is discussed in detail in the next section. First we illus- 
trate in detail the model CMDs made with Salpeter IMF 
{a = 2.35), varying the age range, for the two metallicities: 
Z = 0.004 (Fig. El and Z = 0.008 (Fig.^j). The duration 
of the star formation episode is indicated in each panel as 
the starting-ending time. An ending time of Myr means 
that the star formation is on-going. 

The sharp blue edge of the blue loops is due to the 
theoretical form of the blue loops. In particular, the early 
termination of the 9 M0 model for Z — 0.004 produces 
a bend in the simulated CMDs. This transition is much 
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Fig. 7. Synthetic CMDs for Z ~ 0.008 using the Salpeter IMF (a — 2.35), the appropriate distance, and a constant 
SFR in the interval that is indicated in each panel (beginning-end of the burst). For comparison the observed CMD 
is also displayed on the left. In the lower panels the simulated luminosity functions (solid lines) are compared with 
the observed luminosity functions (dotted lines). In the observed and simulated luminosity functions the color and 
magnitude cuts are applied, and only stars bluer than ([/ — y) = —0.7 and brighter than V — 26 are compared. 



smoother for the Z = 0.008 models. The appearance of the 
CMD suggests that the 9 Mq model for Z = 0.004 should 
have bluer blue edge. The high degree of contamination 
prevents us to draw a robust conclusion, but the U — V vs. 
V CMD seems particularly suited to constrain the shape 
of the loops, and could be applied in stellar systems with 
spectroscopic abundance determination. 

As the star formation lasts longer it is possible to ob- 
serve the excursions of the stars in PMS phases. The red- 
dest parts of the CMD, however, remain very sparsely pop- 
ulated due to the much shorter evolutionary timescales in 
the red with respect to the blue. Less massive (M ^ 15Mq) 
red supergiants evolve less rapidly. Few of them are visible 
in the upper red part of the simulated CMDs in which star 
formation was not terminated. Also, as the star formation 
gets longer, the MS gets populated with stars of smaller 
masses, the total number of which is governed by the IMF. 



In the direct comparison of the simulated and observed 
CMDs, and the corresponding luminosity functions (Fig.El 
and [7)) it can be seen that all the simulations with on- 
going star formation fit the data equivalcntly well irre- 
spectively of the maximum stellar age (Tmax)- The lumi- 
nosity function peak is best fitted with models with longer 
Tmax (^ 30 — 50 Myr), but the differences between various 
models are within the Poissonian errors. Hence, using only 
MS stars {U — V < —0.7) it is not possible to constrain 

On the other hand, the simulations in which the star 
formation has stopped 5 or 10 million years ago, have 
many more stars in the lower MS and around the blue edge 
of the loops than observed. In order to populate the box 
on the MS between magnitudes 23 and 24, which contains 
the youngest MS stars (younger than 14 — 15 Myr) and 
thus the ones with the shortest evolutionary timescales. 
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Table 3. Input parameters and results of the simulated star formation histories of young stellar populations. Columns 
1 to 5 are identifier of the model, slope of the IMF, start and end of the star formation burst in Myr, respectively. 
In the next 8 columns extracted mass, mean star formation rate in the interval Ar, the total number (Ntot) and the 
number of blue plume stars (Nbiue) in the simulated CMDs are given for Z=0.004 and Z=0.008 

Z=0.004 Z=0.008 
model IMF Tmax rmi„ mass (SFR) Ntot Nw„e mass (SFR) Ntot ^biue 
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many more lower mass stars were produced in these sim- 
ulations. Stopping the SF at 2.5 Myr produces a MS lu- 
minosity function in agreement with the observations. 



In the left panels of Fig.|Slwe plot the number of stars 
in the blue plume for all the models with Z = 0.008 (top) 
and Z = 0.004 (bottom). The observed value is shown as a 
green star. Different symbols are used for models with dif- 
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Fig. 8. Left panels: Number of stars in the blue plume (\J — V < —0.7) for V brighter than 25.5 vs. logarithm of star 
formation duration is shown for the models with Z = 0.008 (top) and Z — 0.004 (bottom). The observed number of 
stars within the same limits is indicated with a green star. Dotted lines indicate 3cr Poissonian fluctuation limits in 
the number of observed and simulated stars. The x-axis value for the observed point is arbitrary. Different symbols 
are used for models with different IMF slopes as indicated in the legends. The models for which star formation was 
terminated 2.5, 5 and 10 Myr ago are plotted with open symbols, and the models for which the star formation is still 
on-going have filled symbols. 

Right panels: The ratio of the blue plume vs. total number of simulated stars. The observed point (corrected for the 
foreground contamination) is only a lower limit due to a large uncertainty in foreground contamination. 



ferent IMF slopes as indicated in the legend. The dashed 
lines indicate the 3(t Poissonian deviation from the ob- 
served value to which the random errors of the models 
have been added in quadrature. The models outside of 
the dashed lines in this figure are ruled out. So, the total 
number of MS stars and their luminosity function con- 
strain Train to be at most 2.5 Myr. The impact of the IMF 
on this constraint is discussed in the next section. 

Although PMS stars cannot be used to select the best 
fitting model due to large uncertainties as discussed above, 
in the observed CMD there arc 17 bright (21 < V < 23) 
and blue stars (—0.7 < U — V < most probably in 
the PMS phases. They are found in the area of the CMD 
where there is no contamination by foreground Galactic 



stars. The models that have active star formation only 
over the last 5 — 10 Myr, have one or no stars in this area 
of the CMD. Thus, the star formation in this field must 
have been going on for more than 10 Myr. 

Better limits to the duration of the star formation can 
be obtained by comparison between the number of stars 
in the MS and in the PMS phases. When comparing the 
observed and the simulated number of stars in the red part 
of the CMD, it is necessary to subtract the expected num- 
ber of foreground stars according to the Galactic models. 
The ratio of the blue plume {U — V < —0.7), and the total 
number of stars versus the duration of the star formation 
for the models described in Tab. 01 is presented in the up- 
per right panel of the Fig. |H1 for Z = 0.008, and lower 
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right panel for Z — 0.004 models. This ratio decreases 
with increasing T^ax, because the PMS phases become 
more populated. 

The inflexion point in these diagrams corresponds to 
the maximum look-back time and hence the maximum 
age of the burst Tmax- It is ^100 Myr for the models with 
Z = 0.008 and - 90 Myr for the Z = 0.004 models. Of 
course, this value is the lower limit for the duration of the 
star formation. The older stars are not observed due to 
the incompleteness in the [/-band, hence the ratio of the 
blue vs. total number of stars becomes constant for longer 
star formation durations. 

The observed ratio of the blue vs. total number of stars 
in NGC 5128, which has been corrected for the expected 
amount of foreground contamination, is indicated with a 
star in Fig. |S1 The arrow shows that it is only a lower 
limit due to the large uncertainty in contamination which 
affects in particular the PMS phases. The number of stars 
redder than {U — V) > in the observed CMD is much 
larger than in all the simulations which reproduce cor- 
rectly the number of stars in the blue plume. Since this 
ratio is low, even after the foreground contamination sub- 
traction, it points, independently from the evidence on 
number of MS stars, toward long Tmax- 

3.2. The initial mass function 
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Fig. 10. The comparison of the the VK CMD with stellar 
evolutionary tracks for Z = 0.004 (upper left) and Z = 
0.008 (upper right panel) and masses of 6, 9, 12, 20 and 
40 M0 . In the lower panels the spatial distribution of the 
blue {V - K < 2 and K < 24) stars on the left and red 
{V — K > 5 and 19 < K < 21) stars on the right is shown. 
Note the alignment of the blue stars along the direction 
of the radio jet indicated by the arrow in the lower right 
corner of the left panel. 



The observations of the young stars on the MS can be 
used to constrain the initial mass function for the range 
of masses observed, M ^ 5 Mq. The comparison of the 
simulated and observed luminosity functions along the MS 
is presented for a set of models with on-going SF lasting for 
100, 50, 30, and 10 Myr, as well as for two models in which 
the SF was stopped 5 and 2.5 Myr ago. Fig. [U] compares 
the observed luminosity function (dots) with the models 
for Z=0.008 and four different IMF slopes: Salpeter IMF 
(slope a = 2.35; solid black hue), Scalo IMF (slope a = 2.6 
for M < 3.5 Mq; ScalofTggSl green short dashed hue), IMF 
with slope a = 2 (blue long dashed line), and IMF with 
slope a = 3 (red dotted line). 

The observed slope of the MS between 22.5 <V < 24.5 
is A log N/ AV — 0.47 as derived by least square fitting of 
the observed luminosity function in this magnitude range. 
The results for the models with Z = 0.004 and Z = 0.008 
are equivalent and we plot only Z — 0.008 models in Fig.^ 
The models with IMF slope a = 2.35 provide the best fit of 
the luminosity function for both metallicities, Z = 0.004 
and 0.008. Slightly shallower {a = 2) and steeper (a = 
2.6) slopes for IMF are also acceptable, while the a = 3 
slope produces too many stars in the lower mass bins in all 
the simulations. Note that the models where the SF has 
stopped 2.5 Myr ago are well fitted with shallower IMFs, 
and they would exclude Scalo or steeper IMF slopes. 

These results are confirmed by the tests applied 
to the luminosity functions of the blue stars and are also 
summarized in Fig. |S1 In the left panels the dashed lines 
exclude the models with IMF slope a = 3, because the 



simulated number of blue plume stars is always larger than 
observed (having fixed in the simulations the number of 
stars in a box on the MS). In the right panels the above 
mentioned saturation value (or the inflexion point) of the 
number of PMS vs. the total number of stars depends on 
the IMF slope. The steeper the IMF, the more favoured 
the PMS. Even though the observed ratio would favour 
a = 3.00, this IMF slope is excluded by the excess number 
of stars in the blue plume, which is the more robust con- 
straint. Therefore it seems that the contamination of the 
red part of the CMD is indeed much larger than that given 
by the Besangon counts as already noted by Rejkuba et 
al. pM|l . Additionally, there may be few compact back- 
ground galaxies contaminating the red part of the CMD. 

Our simulations do not include binary stars. The pres- 
ence of binaries would increase the IMF slope by an 
amount that depends on the slope and the fraction of stel- 
lar mass in binary or multiple stellar systems. According 
to Sagar & Richtler |1991) if the fraction of binaries is 
more than 50% the IMF slope increases by 0.3 to 0.4. As 
the binary fraction in our field is not known, it is possi- 
ble that the actual IMF slope is somewhat steeper than 
a = 2.35. 

4. Young stellar population in the VK CMD 

Fig. [TUl shows VK optical-near IR CMD for this field con- 
structed from the VLT FORSl V^-band and ISAAC Ks- 
band photometry published by Rejkuba et al. (2001). Most 
of the stars are distributed along the red giant branch and 
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Fig. 9. The observed MS luminosity function (cyan dots) is compared with synthetic luminosity functions for different 
SFHs for Z = 0.008 and four different IMF slopes (a = 2.0, 2.35, 2.6, and 3.0) as indicated in the legend. Synthetic 
luminosity functions correspond to models calculated adopting the appropriate distance, reddening and a constant SFR 
in the interval that is indicated in each panel (beginning-end of the burst). In the observed and simulated luminosity 
functions color and magnitude cuts are applied, and only stars bluer than U — V — —0.7 and brighter than V = 26 
are compared. 



are thus older than ~ 1 Gyr. There is also a large num- 
ber of bright AGB stars, including several hundred Mira 
variables, located above the tip of the RGB (Rejkuba et 
al. nnn^ . some ~ 110 stars with colors 2 <, V - K 4: 
are expected to belong to foreground Galactic population. 
Several dozens of stars in the bluest part of the diagram 
{V — K < 2) do not appear in the control halo field (see 
Rejkuba et al. I2001|l . They have very small photometric 
errors (90% has a < 0.2), and are bright in the optical. All 
but two stars are detected in the ?7-band images as well. 
We plot the spatial distribution of the blue {V — K < 2 
and K < 24) and red {V - K > 5 and 19 < K < 21) 
stars in the shell field in Fig. ^1 While the reddest ob- 
jects spread evenly over the field, the majority of the blue 
stars are aligned in the direction of the jet just like in the 
optical CMDs. 



The low number of young stars observed in the much 
smaller field of view of ISAAC prevents us from drawing 
conclusions from their observations. The comparison with 
the Padova stellar evolutionary tracks (in the upper pan- 
els of Fig. nni) for two different metallicities Z = 0.004 
(left) and Z = 0.008 (right), indicates again a very small 
difference between the models. 



The VK CMD gives complementary information 
about the older stellar population in the halo. However, 
due to large incompleteness in the red parts of this CMD 
(see Rejkuba et al. I2001|l . the SFH of the old stars in 
NGC 5128 will be investigated with deeper HST+ACS 
observations. 
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5. Discussion and Conclusions 

The recent star formation history in the north-eastern halo 
field in the giant elliptical galaxy NGC 5f 28 was investi- 
gated with the means of synthetic CMDs. The compari- 
son of the isochrones and the data has to be done with 
great care, using the correct filter transmission curves to 
compute a consistent set of bolometric corrections. The U- 
band is known to be the most difficult to reproduce due 
to the presence of the atmospheric cut-off, to the pres- 
ence of the Balmer break within that band, and to strong 
blanketing by metals. Additionally due to large differences 
between the tJ-band filters and a possible miss-match in 
stellar properties (logg,re// and [M/H]) between the cal- 
ibrators and program objects, it is possible to introduce 
systematic errors if the calibrators do not rely on the 
same photometric system as the observations. This affects 
mostly our reddening estimate and very little other re- 
sults. 

The new reddening estimated by matching the locus of 
the observed MS in the U — V plane with the models, us- 
ing the corrected bolometric corrections for U passband 
and assuming the Galactic extinction law (Cardelli et 
al. ll989ll is E(i3 -V) = 0.35, significantly larger than pre- 
viously obtained E(B -V)^ 0.15 (Rejkuba et al. [2001,1 . 
It should be noted, though, that Rejkuba et al. estimate 
is obtained from the comparison of mostly foreground 
stars with isochrones, and is in fact in good agreement 
with E(B -V) ^ 0.11 from Schlegel et al. maps. 
This implies a rather significant internal extinction in this 
field. While some internal extinction in the areas of re- 
cent star formation could be expected, we would like to 
stress that this result has large uncertainty due to above 
mentioned difficulties to reproduce correctly the J7-band. 
The maximum possible systematic error introduced by 
the calibration to the Landolt photometric system is of 
the order of 0.04-0.08 in the [/-band for the bluest stars. 
This translates to a possible over-estimate of reddening of 
E(i3 — V) = 0.07 — 0.14. Assuming the maximum possi- 
ble systematic error, the internal extinction in this field is 
reduced to about E(B ~V) ^ 0.10. 

Unfortunately, on the basis of the U - and y-band pho- 
tometry alone it is not possible to constrain the metallicity 
of the stars due to reddening-distance-metallicity degen- 
eracy of the models. However, the comparison of the syn- 
thetic and observed CMDs gave a number of quantitative 
constraints to the recent SFH in the studied halo field. 

We have constructed simulated CMDs of young stars 
in the halo, having fixed the number of stars in a box on 
the bright portion of the MS, for two metallicities Z=0.004 
and Z=0.008. The parameters that were varied were r^ax 
(the starting epoch of the SF), Tmin (the ending epoch of 
the SF), and a (the slope of the IMF). 

The blue MS stars and their luminosity function con- 
strain Tmin to be smaller than ^2.5 Myr and exclude the 
IMF slopes steeper than a 2.6. The ratio of the blue 
vs. total number of stars favours steeper IMF slope and 
longer star formation duration {T^ax ^ 100 Myr). This is 



a robust result and is not affected by the systematic errors 
in calibration, as the slope of the MS does not change. The 
exact ratio of the blue vs. total number of stars may be 
affected by the systematics, but not so much to change 
our conclusions. 

From the simulations with Tmin = and Tmax = 
100 Myr, the total extracted mass in stars in the range of 
masses 0.6 < M < 100 M© is ~0.4, 0.7, or 1.3 x 10^ Mq, 
for IMF slopes of 2 and 2.35, or 2.6, respectively. The SFR 
associated with these IMFs and mass limits are approxi- 
mately 0.004, 0.007, and 0.013 Mq yr"\ for the three IMF 
slopes in the above order, and it is independent of Tmax- 
Moreover the smooth appearance of the MS, the absence 
of gaps and breaks in the slope of the MS luminosity func- 
tion, indicates that most probably no major fluctuations 
in the SFR have occurred over the sampled time interval. 
The reported values of the SFR depend on a combination 
of the MS and PMS stars, for which the possible system- 
atic errors in magnitudes (that go in different direction for 
these two groups) partially cancel out. We estimate that 
our results are not changed strongly by the systematics in 
calibration. 

At the distance of 3.6 Mpc, the observed field surveys 
a projected area of approximately 45 kpc^, so that the 
derived star formation rates range from ^ 9 x 10^^ to 
~ 3 X 10~* M.Q yr~^ kpc~^. This is lower than the typ- 
ical SFR quoted for the solar neighbourhood (4.2 x 10"^ 
Mq yr~^ kpc~^; Timmes et al. I1995|l . and is on the low 
side of the SFRs derived via theoretical simulations for 
dwarf irregular galaxies in the Local Group (ranging from 
approximately 10~^ to 10~^ Mq yr~^ kpc~^; Marconi et 
alUnnSl TolstoyCnni Callart et aL ll996a|l . For example, 
the star formation rate in WLM dirr within the last mil- 
lion years was 1.3 x 10"^ Mq yr'^ kpc"^ (Dolphin OUnUjl . 

On a more speculative side, it is possible to estimate 
the maximum duration of the future star formation. The 
fuel for this recent star formation is provided by the large 
HI gas cloud in the vicinity of the field (Schiminovich et 
al. 11994(1 . Assuming the constant SFR and a 100% effi- 
ciency in converting gas to stars, we can estimate the gas 
depletion timescale as: 



In the north-eastern shell very close to where the star for- 
mation is observed, 2.1 x 10^ Mq of HI have been detected 
(Schiminovich et al. I1994 II as well as 1.7 x lO'' Mq of H2 
(Charmandaris et al. |2)00). At the constant star forma- 
tion level of 0.006 Mq yr-^ it would take ~100 Myr to 
use up all the molecular gas and ^^230 Myr to use all the 
hydrogen. 

Our results place a lower limit of ~ 100 Myr to the 
last merging event and a lower limit to the time of the 
formation of the jet in Cen A. To set better constraints 
on the starting epoch of the last SF episode deeper CMDs 
and a more robust determination of the foreground con- 
tamination are necessary. 
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Appendix A: Random extraction of a star with a 
given mass and age 

We use a Monte Carlo algorithm to randomly select a mass 
and an age according to a given distribution. The mass is 
extracted according to a single slope (e.g. Salpeter IMF; 
SalpeterCnnni), 

dN = AM-^dM (A.l) 

or double slope initial mass function (e.g. Gould IMF; 
Gould et al. 11997(1 with slope ai within a mass range 
Ml < M < M2 and ^2 within M2 < M < M3. Here 
"A" is the normalisation factor that depends on the up- 
per and lower mass limits used in the models. The choice 
of the slope is a free parameter. 

The age distribution is adopted to be either flat, rep- 
resenting a constant star formation rate, 

*(t) C for T„i„ < t < Tmax (A.2) 

or exponentially decreasing, 

^{t) = Ce-*/" for r™„ < t < r^ax (A.3) 

with given limits for the age of the oldest stars {Tmax) 
and for the duration of the star formation episode (Ar = 

We have used a single metallicity simulator to produce 
synthetic CMDs for the young stellar population. This is 
a good approximation, because the spread in ages of the 
young stars (few x 10 Myr) is too short for a significant 
metal enrichment to occur. However in the simulation of 
the optical-near IR CMDs that are dominated by the old 
and intermediate-age population, evolution of the metal- 
licity would have to be assumed. 

Appendix B: Interpolation of the stellar 
evolutionary models 

We have adopted two different sets of models: Padova 
models (Fagotto et al. 11994;) and models computed 
by Cassisi & Salaris H1997fl and extended by Bono et 
al. H1997al 119976 1. For the simulations of single metal- 
licity star formation histories used to compare with the 
youngest stars in NGC 5128, we use the Padova models 
by Fagotto et al. l(1994|l . These models have tracks for 
stellar masses ranging from 0.6 to 120 Mq. The newer set 
of Padova models (Girardi et al. I2000fl are not computed 
for stars with masses larger than 7 Mq such as those that 
are observed in the north-eastern halo field of NGC 5128. 

Given the extracted mass and age, the luminosity and 
temperature of the star were determined via the linear 
interpolation in a set of evolutionary tracks. Each track 



consists of a number of points (age, mass, log(L/L0), 
log(Toff)). We have divided each track into portions ac- 
cording to stellar evolutionary phase, and interpolated lin- 
early between the adjacent tracks within the same portion. 
Typically, the portions of the tracks are: 

— Main sequence (MS) - ending either at core H exhaus- 
tion or at the start of the first runaway expansion 

— Subgiant branch - ending at the base of the RGB 

— RGB - ending at the tip of the RGB 

— Horizontal branch - ending either at the core He ex- 
haustion or at the base of the AGB 

— Early AGB - ending at the first thermal pulse 

For the high mass stars (M > 20 Mq), whose evolution 
is affected by mass loss, we divide the evolution in two 
major phases (MS and post main sequence), and choose 
conveniently the description of the post main sequence 
(PMS) phase, inserting further sub-phases when needed. 

In most of the cases the log(L) and log(Toff) of the 
extracted star are interpolated between the values of ad- 
jacent masses read at the same fractional age in the ap- 
propriate evolutionary phase: 

tph = {t- to)/{ti - to) (B.l) 

where to and ti are respectively the ages at the beginning 
and at the end of the appropriate evolutionary phase and 
t is the age of the extracted object. In few evolutionary 
phases when luminosity and/or temperature presented a 
maximum or minimum at variable tph , we inserted a break 
tph{tph,br) and effectively divided that phase into two sub- 
phases. Log(L) and log(Toff) of the extracted star at the 
extracted age are then found by interpolation at: 

tph/tph.br for ages tph < tph^tr (B-2) 
{tph — tph,br)/{^ — tph,br) for ages tph > tph,br (B-3) 

When the interpolation has to be done also in metal- 
licity, the simple linear interpolation algorithm described 
above is not enough. It is necessary to perform a bilinear 
interpolation, because the stars of the same masses, but 
different metallicities have different ages. Thus for each 
extracted mass one has to first interpolate in log(L) and 
log(Toff ) for a given age, and then again for a given metal- 
licity. In order to simplify the problem Zoccali et al. H2003|l 
adopted a set of isochrones as a starting point. Using the 
tabulated values of log(L) and log(Tcff ) for stars of given 
mass in an isochrone the interpolation is done linearly 
within isochrones of the same age and different metallic- 
ities. The Padova models are not adequate for these sim- 
ulations due to the very low number of points in each 
isochrone (approx. 30-55 points for each track) and it is 
better to use the Cassisi models (Cassisi & Salaris [T997\ 
Bono et al. ll997alll997b|l which have much more detailed 
evolution with more than 600 points per track. Given a 
mass and metallicity, the luminosity and temperature of 
the star are determined through the linear interpolation 
in a set of isochrones. The interpolation in log(L) and 
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log(Teff) is made for the mass that has the same evolu- Kurucz, R.L., 1993, lAU Symposium 149: "The Stellar 
tionary phase of the extracted mass: Populations of Galaxies", B. Barbuy, A. Renzini, eds., 

Dordrecht, Kluwcr, p. 225 
ph = m/mi for objects on the main sequence (B.4) Landolt, A.U., 1992, AJ, 104, 340 
ph = mo + (m- mo) /(mi - mo) for later phases (B.5) Marconi, G., Tosi, M., Greggio, L. & Focardi, R, 1995, AJ, 

109, 173 

where mo and mi are the masses at the beginning and the Montegriffo, P., Ferraro, F.R., OrigUa, L., Fusi Pecci, F., 1998, 



end of the cvohitionary phase. The evolutionary phases 
are as before: main sequence, subgiant branch, RGB, core- 
helium burning and AGE. 
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